Topology in lattice QCD 



A. Di Giacomo 

. Dipartimento di Fisica and INFN, 

Q> [ 2 Piazza TorriccUi 56100 Pisa, Italy 

digiacomopi.infn.it 

> 

o 

QQ ■ Abstract: The status of topology on the lattice is reviewed. Recent results show that the topological 

susceptibility x can be unambigously determined. Different methods, if properly implemented, give results 
consistent with each other. For SU{3) the Witten - Veneziano prediction is confirmed. Preliminary results 
for full QCD are presented. The problem there is that the usual hybrid montecarlo algorithm has severe 
difficulty to thermalize topology. Possible ways out are under study. 
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1 Introduction 

Topology plays a fundamental role in QCD. The key equation is the anomaly of the Ua{^) axial current ||l| 



d^j'Jx) = 2NfQix) (1) 



J^^ = X]/2:i '^fl^ltii'f is the singlet axial current. Nf is the number of light flavours, and 



Q(^) = ^ E G^G"'"' (2) 
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' the topological charge density. = ^e^i/po-G"'"^ is the dual tensor field strength. At the classical level 

d^J^ = 0: the right hand side of eq.(|^) comes from quantum effects, whence the name of anomaly. Q{x) 

is related to the Chern current i^^(x) = if^e'^^'^'^A^ (dpA'i^ - f as follows 

a^X^ = Q{x) (3) 

Eq.(|l|) can then be written 

d'^{4~2NfK^)^0 (4) 

from which Ward identities are derived^. 

A number of physical consequences follow from eq.'s (|^), (^). 
1) The Ua{^) problem. This problem goes back to Gellmann's free quark model, from which the sym- 
metries of hadron physics were abstracted. In that model UA{i) is a symmetry, i.e. 9^J^ = 0. This 
implies either the existence of parity doublets in the hadron spectrum, if the symmetry is Wigner, or, if 
it is Goldstone, Weinberg's inequality!^ for the 77' mass to„' < y/Smyr- Neither is true in nature. 

In QCD C/a(1) is not a symmetry, due to the anomaly (|^) and in principle there is no Ua{^) problem. 
The anomaly can however explain the high value of the ry' mass, as suggested by a 1/Nc expansion. The 
idea is that in the limit A'c — > 00, with g^Nc fixed the theory preserves all the main physical features of 
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QCD, like confinement Q. Non leading terms act as a perturbation on this limiting model. As Nc — s- oo 
the anomaly disappears, being 0{g'^) [eq.(2)], and the j]' is the Goldstone boson of the spontaneously 
broken [7^(1)- The effect of the anomaly is to shift the Goldstone pole from zero mass. This shift can be 
computed by use of the Ward identities (|^) giving[Q 

—fx = K + mi^' - 2tok (5) 
X = x(9^)|^_Q is the topological susceptibility of the vacuum at A^^ = 00, or 

X{q^)= I d^xe"^^0\T[Q{x)Q{0))\Q) quenched (6) 



Nc = 00 implies quenced approximation, fcrmion loops being 0{g^Nf). 
Eq.(|) predicts, within 0{1/Nc) accuracy 

X= (ISOMeV)'^ (7) 

To be definite a prescription must be given for the singularity in the product Q{x)Q{0) in eq.(||) as x — > 0. 
The prescription which brings to eq.(||)is|^ |j 

X = J d\x - y)dldyMT{K,{x)K,{y)m (8) 

By the prescription (|^) 5 like singularities in the product K^{x)K^{y) as x y are eliminated after 
integration. In any regularization scheme only the multiplicative renormalization of Kf^ is thus left. 
Eq.(^ can be verified on the lattice, where Xquenched can be computed. This computation provides a 
check of the expansion, which is a fundamental issue. 

2) The behaviour of x across the deconfining phase transition is relevant to understanding the structure 
of QCD vacuum|]. 

3) The Ua{^) Ward identities predict that in full QCD x behaves in the chiral limit linearly in the quark 
mass ^ 

X-^C^'^fi'f^f) (9) 
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This determination in 



4) Another important quantity to determine on the lattice is x' = df^x(9^) 

full QCD is relevant to understand the spin content of the proton Q. 

5) The spin content of the proton can be determined by use of Eq.(l), as will be discussed in detail in 
the following. 

A lattice regulator Ql{x) of the operator Q{x), eq.(||) is needed. The regularized matrix elements 
are then determined by numerical simulations, and out of them the continuum physical quantities can be 
extracted by proper renormalization. 



2 Defining Q{x) on the lattice. 



According to the general rule any gauge invariant operator on the lattice Ql{x), such that in the formal 
limit a — !■ 

QL{x) = Q{x)a^ + Oia^) (10) 

is a possible regulator of topological charge density. 

A large arbitrariness by higher order terms exists in the choice of Ql , which can be used to improve 
the operator in the sense which will be discussed below. A prototype choice for Ql is 
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(11) 
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where 11'^^^ is the usual plaquette operator in the plane Q^^^ obeys eq.(|To|), and differs by 0{a^) 

from any other choice. We will make use of a recursive improving of the operator (pl])[p[. We shall define 
q''1\x) (i = 1, 2 . . .) by a formula similar to (pi]), with n|'J'^n^°j replaced by Il'^^^Il'^^y For i — I each link 
is replaced by a smoothed link, for z = 2 each link of 11(1) is smoothed again and so on|^. An alternative 
definition is the so called geometrical charge]^ |l^, which again obeys the constraint eq.(|lO|). 

Any Ql, in the limit a —> 0, will be a mixing of all the continuum operators with the same quantum 
numbers and with lower or equal dimension. In quenched QCD the only pseudoscalar with dimension 
< 4 is Q{x) itself and therefore Q 

Ql{x) ~ ZQ{x) (12) 

Z can be determined by a non perturbative procedure known as "heating" ^] . The idea is to measure 
Ql on a state with a definite value of Q, which then determines Z. To do that an istanton is put by 
hand on the lattice, and quantum fluctuations at a given value oi (3 = 2N/ are produced numerically 
by the usual updating algorithm. A plateau after a few heating steps signals thermalization of these 
fluctuations, while the topological content of the conflguration takes a much longer time to be changed, 
as can be directly tested: the values of Z for different operators at different /3's are shown in Table 1. 





Op 


XL ■ 10-5 


M ■ 10-5 


M/XL 


Z 


X^/'' MeV 




Qo 


2.72(6) 


2.50(15) 


.9 


.12(4) 


167(36) 


13 = 5.9 


Ql 


2.48(5) 


.88(6) 


.35 


.36(2) 


178(6) 


= 1.74(4) GeV 


Q2 


3.51(7) 


.71(6) 


.2 


.48(2) 


175(6) 


Nc^3 


Qo 


2.14(4) 


1.97(10) 


.9 


.18(4) 


196(40) 


13 = 6.1 


Ql 


1.12(2) 


.47(3) 


.42 


.41(2) 


178(6) 


= 2.17(5) GeV 


Q2 


1.39(2) 


.33(2) 


.24 


.54(2) 


175(6) 


Nc = 2 


Qo 


3.38(8) 


2.93(7) 


.87 


.18(6) 


208(23) 


13 = 2.44 


Ql 


3.26(8) 


1.45(4) 


.44 


.405(13) 


197(8) 


a-^ = 1.92 GeV 


Q2 


4.92(14) 


1.34(5) 


.27 


.568(16) 


197(8) 




Qceo 


69.7 ± 1.1 


52(2) 


.75 


.918(57) 


231(13) 




Qo 


2.32(5) 


2.20(3) 


.95 


.240(26) 


197(27) 


13 = 2.57 


Ql 


1.01(5) 


.44(2) 


.45 


.507(9) 


200(8) 


= 2.91(9) GeV 


Q2 


1.16(6) 


.117(5) 


.16 


.675(8) 


198(7) 




Qceo 


16.6(3) 


13.26(23) 


.80 


.937(26) 


228(12) 



Table 1. xl, M,x, Z,a for quenched SU{3) and SU{2) at various (3. 



3 Measuring x 

The lattice susceptibility is deflned as 

X 

XL is a positive quantity by definition. However in the euclidean region, for x ^ 0, {Q{x)Q(0)) is a 
negative quantity, by refiection positivity, Q being odd under time refiection. This can be checked on 
the lattice and the result is shown in fig.l. If the operator is smeared over a region of size s, one expects 
{Q{x)Q{0)) to be negative at distances \x\ > s. The peak at x = is essential to make x positive, and 
hence the prescription for the product at x = is essential. In general, as a consequence of operator 
product expansion around x = 0, xl will mix with the continuum operators with the same quantum 
numbers and lower or equal dimension[^3| 

XL = ZV + M(/3)+0(a') (14) 
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with 



M{(3) = Z(/3)a^ 



//3(5) 
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Z{f3)a'^{mip{x)il;{x)) + P{/3){I) 



(15) 



In the quenched approximation the second term in eq.(|l5|) will be absent. 

The first term in cq.(p^ corresponds to the prescription and renormalizes multiplicatively. Since 
/ d^xd^K^^ = Q the prescription (^) implies that x niust be zero on the sector Q — 0. M{(3) can 
be obtained by measuring xl on that sector. This is done by the same heating technique [T^ used to 
determine Z: a sample of configurations belonging to the Q = sector is produced by heating the flat 
(zero field) configuration, in such a way that the topological charge is not changed. From eq.(Qj) then 



XL - M{(i) 
a\(3)Z^{(3) 



(16) 



All the quantities in the r.h.s. of eq.(p^ depend on the choice of the operator, and/or on the choice of 
the action, but the result must be independent of them. This appears from Table 1, and from fig. 2, 3. 

A good operator is such that |Af(/3)|/|xL| ^ 1 so that most of the observed signal is physcal. Z ~ 1 
is also desirable. The quality of different Ql's can be appreciated from Table 1. 

The behaviour of x at deconfinement is shown in fig. 2 iQ, where SU{2) and SU{3) can be directly 
compared. 

The existence on lattice of a nontrivial Z was first realized in ref.13. Eq.([l6|) was first introduced in 
ref.ll, where, however, the renormalization constants were determined by perturbation theory. 
The heating technique jl^ finally allowed a non perturbative determination of them. 

Fig. 1 Fig. 2 





fig.l Correlation {Q{x)Q{0)) for the geometric 
topological charge density. 



Fig. 3 



fig. 2 X for quenched SU{3). Diamonds, circles and 
squares corrspond to 0,1 and 2 - smeared operators. 

Fig. 4 



OSU(2)dala 
I □ SU(3) dala | 



T/T. 



fig. 3 X for quenched SU{2). 



fig. 4 X across the deconfinement transition for SU{2) 
(circles) and SU{3) (squares). 
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4 Full QCDi^. 



We have used the same technique to determine x in fuU QCD with staggered fermions at mg — 0.01 and 
TOQ = 0.02 at /? = 5.35. The lattice was 16^ x 24. At this (3 value a = 0.115(2) fm. 

We do not have enough precision to test the dependence on mq [eq.(^)]. At am — 0.01 we get 
X = (110 ± 6 ± 2)"* MeV'' The first error is obtained by jake-knife technique. The second comes from the 
error in the determination of the scale. The expected value is x — (109 MeV)^. 

As we will soon discuss these errors are higly underestimated. 

A similar analysis can be done for x' ■ On general grounds 

= Z^x'a^ + M'ifi) + 0{a') M'{(3) = P\(3){I) (17) 

x' only mixes to the identity operator, since no other gauge invariant operator of dimension < 2 exists. 

The technique used to determine x' is to improve the operator enough, so that M'{(3) will be neglige- 
able as compared to the first term. Z is known from the analysis of x- 

A preliminary result is 

X' = (258 ± 100) MeV^ Vx' = (16 ± 3) MeV (18) 

to be compared with the value computed from SVZ r ules§ ^ (22.3 ±4.8) MeV Here the error coming 
from the normalization is much smaller than the statistical error. 

Again here the error is underestimated. The reason can be in the history of the topological charge 
along the updating process is shown. The hybrid montecarlo algorithm is very slow in thermalizing 
topology [p^. The sampling is very bad and corresponds in fact to a much smaller number of independent 
configurations. 

The same incovenience affects the determination of the spin content of the proton which will be 
discussed in the next section. The real error is then larger than the one estimated on the present 
ensemble of configurations. 



5 The spin content of the proton. 

The matrix element of the singlet axial current between proton states can be written as 

{0s'\Jlmps) - u{0s') {Gi(fc2)7'7M + G2{e)l'k^) u{ps) (19) 

k = p — p' is the momentum transfer. Gi(0) = AT, is related to the integral of the spin dependent 
structure function gi(x,q^) of deep inelastic scattering of leptons off protons. In naive parton picture 
AE = Au + Ad + As is the fraction of theproton spin carried by the quarks and is expected to be 
AE = 0.7. Experiments [p^ show that it is much smaller, AE = 0.21 ± 0.10, and this fact is usually 
referred as spin crisis. By use of the anomaly equation 

A{k'^) = Gi(/c^) + (fc^/mjv)G2(fc^) can be extracted from a measurement of the matrix element in the 
left hand side of eq. ([l9|) . 

In quenched approximation G2{k^) has a pole at k'^ = 0[|9j. With dynamical quarks this is ot the 
case and, for sufficiently small k'^ Gi(0) or AS[^ can be determined. 

Our preliminary result is AE = 0.05 ±0.05 Again the error is underestimated due to the bad sampling 
of topological charge discussed above. 
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6 Discussion. 

Our conclusions are the following. 

A reliable determination of the topological susceptibility x is possible on the lattice. All existing 
definitions of lattice topological charge density give results consistent with each other if proper renor- 
malizations are performed. The value of x for SU (2) is larger than the prediction!^ but the value for 
SU{?>) is consistent with it. The value of m^/ is well explained by the anomaly. Above the deconfining 
transition x drops to zero, more rapidly for SU{'i) than for S't/(2)|lJ]. In full QCD preliminary results 
are in agreement with expectations for x ^^'^ x' ^-i^d for the spin content of the proton. However a better 
sampling of topological sectors is needed. In fact the incapability of the hybrid montecarlo to thermal- 
ize topology could affect not only the measurement of quantities directly related to it, like the one we 
considered, but any other measurement on the lattice. Work is in progress to overcome this difficulty. 

I am grateful to my collaborators B. AUes, G. Boyd, M. D'Elia, E. Meggiolaro, H. Panagopoulos, with 
whom most of the results presented where obtained. 
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